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Abstract. Isospin-breaking effects in the baryonic sector are studied in the framework of a medium-modified 
Skyrme model. The neutron-proton mass difference in infinite, asymmetric nuclear matter is discussed. In 
order to describe the influence of the nuclear environment on the skyrmions, we include energy-dependent 
charged and neutral pion optical potentials in the s- and p-wave channels. The present approach predicts 
that the neutron-proton mass difference is mainly dictated by its strong part and that it strongly decreases 
in neutron matter. 
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1 Introduction 

The medium dependence of isospin-breaking effects be- 
longs to one of the fundamental questions in nuclear physics 
[IK^lIS]. In particular, the neutron-proton mass difference 
in nuclear matter Z\rn*p is an interesting topic of nuclear 
astrophysics relevant to the evolution of the universe at an 
early stage [31[S]. Furthermore, it is also important for the 
description of the properties of mirror nuclei [6] , the stabil- 
ity of drip-line nuclei [7] and the transport in neutron-rich 
matter induced by heavy- ion collisions [8J . Although there 
exist various publications dealing directly with the density 
dependence of the neutron-proton mass difference [MTUl 
fTTlfTlfT3l[T4l[T5lfT6lfT7l and its impHcations for asymmetric 
nuclear matter and finite nuclei properties P^finil^Ull^ 
[^^3, 24, 25, 26, 27, 28, 29, 30, 31, 32J. this quantity is stiU 
not well understood. Quantitatively and even qualitatively 
the predictions about the behavior of the neutron-proton 
mass difference in nuclear matter change from model to 
model. 

Skyrme-soliton models have the inherent advantage 
compared with other hadronic models that they are based 
on chiral input and that they treat the properties and 
interactions of the nucleons on an equal footing 33,34 . 
In this context we recently studied isospin-breaking ef- 
fects in the baryonic sector of a medium-modified Skyrme 
model |35| . by focusing on the single hadron properties 
in the nuclear environment rather than on the properties 
of the system as a whole. The approach predicted that 
the neutron-proton mass difference changes in an isospin- 
symmctric nuclear environment by a very small amount. 



The isospin-breaking leading to this result was only due to 
a modification of the mesonic sector of the Skyrme model, 
originally introduced to generate the strong neutron-proton 
mass splitting in free space (in addition to the electro- 
magnetic one) [36J. However, when the nucleons are em- 
bedded in an isospin-asymmetric environment, additional 
medium effects can be expected. To evaluate the latter in 
the present work we will consider the nucleon properties 
in homogeneous, infinite and isospin-asymmetric nuclear 
matter. The Lagrangian of the present study is a gen- 
eralization of the in-medium Skyrme-type Lagrangian of 
Ref . [35] . In addition to the strong isospin-breaking in the 
mesonic sector it explicitly takes into account the differ- 
ent influences of the isospin-asymmetric environment on 
the charged (tt^) and neutral (tt") pion fields via a built- 
in energy dependence in the s-wave pion- nucleon (tt^N) 
scattering lengths [37,38,39,40,41 and the p-wave opti- 
cal potential [32] • This additional energy dependence also 
alters the predictions for iso-symmetric matter consider- 
ably. 

The paper is organized as follows: In Sect. [2] we for- 
mulate the model including the optical potential input 
and medium modifications. Section [3] discusses the clas- 
sical Lagrangian and the pertinent equation of motion. 
In Sect, m we present the quantization procedure and the 
final expressions for the the strong and electromagnetic 
part of the in-medium neutron-proton mass difference. 
The results of the calculation are reported and discussed 
in Sect. [51 In Sect. [Hi our conclusions are summarized and 
an outlook to future studies is given. Moreover, for clar- 
ification, two appendices are added. In the first one, we 
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discuss in detail the setup of the classical equation of mo- 
tion in the presence of isospin-breaking terms. The second 
appendix is devoted to peculiarities in the construction of 
the charges and magnetic moments. 



2 Formulation of the problem 

2.1 Medium modification of the model 

As in our previous work [35j . we start with a generalized 
Skyrme-type Lagrangian which incorporates an explicit 
isospin-breaking term in the mesonic sector: 



XSB 



'C2 = -^Tr(i^L^) , 



-gxSB 



1 



32& 

i?2 



2 Tr [Z/^, , 



(1) 
(2) 

(3) 
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{Tr[(C/-l)Xi(C/t-l)] 
-Tr[(C/-l)r3X^(C/t-l)r3]}, (4) 



where Einstein's summation convention is always assumed 
(if not specified otherwise). = U^df^U is given in terms 
of the chiral SU{2) matrix U = exp(2zTa7ra/F^), where 
ttq (a = 1,2,3) are the Cartesian isospin-components of 
the pion field. F^^ = 2/7r is the pion-decay constant, while 
e is the dimensionless Skyrme constant. Finally, M± = 



(m^-t ± m^o)/2 is defined in terms of the masses of the 

charged and neutral pions. As in Ref. [3S] we insist on re- 
producing the empirical (isospin-averaged) masses of the 
nucleon and delta, tun = 938 MeV and Ma = 1232 MeV, 
in free space (density p — 0) and without isospin break- 
ing term (A^_ — 0). Furthermore, as input for the free 
mass of the neutral pion we take the PDG- value [43] : 
mTrO = 134.977 MeV. These choices induce the values 
= 108.11 MeV and e = 4.826. Following Ref. [35| the 
mass of the charged pions m^± is extracted as a variational 
parameter m^± = 135.015 MeV from the fit to the empiri- 
cal value Am^p^"^ =1.29 MeV of the neutron-proton mass 
splitting in free space. Note that the dominant electromag- 
netic contribution to ■m^± — m^o is beyond the scope of 
the model. 

The generalized pion mass term Q, which was origi- 
nally proposed by Rathske [36] , can be rewritten as 

Cs^SB = -§{Tr [{U ~ l)mUU^ - 1)] 

+ J2MraU)MlTT{TaU^)\ , (5) 

a=l ' 

which is convenient for our up-coming modifications. When 
the pion fields are small. 



U = exp 



2iT ■ TV 

F^ 



2iT ■ TV 



(6) 



the Lagrangian ([T|) reduces to the Lagrangian for free pion 
fields 



+ i(9,^°9'^7r°-^V2„^0) , 



V2 



(7riT«7r2), 7r"=7r3. 



(7) 
(8) 



In the medium the analog of the Lagrangian ([7]) reads 





A=±,0 
£low 



(9) 



where /Z" and All are linear combinations of the self en- 
ergies of the charged pions: 



n° = -(n- 

2^ 

An = l{n- 



(10) 



As in Ref. [35j the medium-modified version (always mar- 
ked by an asterix) of the Lagrangian ([1]), can be defined 
as 

C* =C2+U + Cl^SB. (11) 
where for the general case of asymmetric matter £g^gB is 
given by formulaQ 



^gxSB = -^{Tr[(C/- 1)K + 77°)(C/t - 1) 



E Tr(raC/) 



a,b=l 



zeafc3^^/2] Tr(TbC/^)} 



(12) 



It is easy to check that the Lagrangian pT|) reduces to 
the Lagrangian ([9]) for the expansion ([6]) as well as to 
the medium-modified Lagrangian of Ref. [351 for the case 
of isospin-symmetric matter, 77+ ~ TI^ — (in the 
parameterization of Ref. [35j). 



2.2 Parameterization of the optical potentials 

The polarization operators of the charged pions can be ex- 
pressed in terms of energy-dependent pion-nucleus optical 
potentials [32] as follows: 



11 f (c^, r) = -4^6± {uj, r) =xt{^,r) 

47rc^(w, r) 47rti; 



(13) 



i7±(c^,r) = ^ ' , • V - (vV(w,r)) 

P ' ' l + 47r5r'c±(w,r) 2m n ^ ' 

^ From now on m„ stands for the mass of the neutral pion, 
i.e. mir = m^o. 
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6±(c.,r) ^ (6S«Hp(r)T6iM<5p(r))r?, 



P{r) = pn{r)+ Pp{r) 
Sp{r) ^ p„{r) ~ pp{r) 



(14) 
(15) 

(16) 

(17) 
(18) 
(19) 



where p„ and pp are the neutron and proton densities, 
respectively. Note that additional V^p and 'V^Sp terms 
are included in the p-wave optical potential since they are 
needed for the description of realistic pion-nucleus scat- 
tering data [42] , 

The chiral expansion of the off-shell pion-nucleon scat- 
tering amplitudes at vanishing pion three-momentum leads 
to energy-dependent s-wave isoscalar and isovector scat- 
tering lengths, 6o('^) a-nd &i(w), respectively. The quanti- 
ties cq{lo) and ci(aj) are the corresponding p-wave scat- 
tering volumes, whereas &q^(w) is the effective isoscalar 
scattering length (see Eq. ([22]) ). The correlation parame- 
ter g' , which renormalizes the pion dipole susceptibility, is 
fixed at g' = 0.47. 

Within the counting scheme of pion-nucleon chiral per- 
turbation theory from Refs. [571155] and based on input 
from these references &o(^) and hi(w) can be expressed at 
order 0{ml) as [MlllQ] 



6o(w) 



1.206to"1 



bp , 
Attt] ^ 

1 LU 



1 ~ m-^LO^] 



1 - m-^LU^) 



1 + 



(20) 



7CJ 



ph 



ph. 



4^'? 2/2 

_ 1.115m-i + 0.143m;3,,3) 

Anrj 



(21) 



Here a^N = -^Ciml - 9g^m^/647r/2 « 45 MeV is the 
pion-nucleon sigma term, whereas the other parameters 
correspond to the "range term" (3 = g\ /4m — 

2C2 — 2C3 w 0.541m~^ and also 7 = (gA7r/7r,ph/wN)^ + 
ln(2ylc/m7r) ~ 2.523. The axial-vector coupling constant 
gA = 1.27 and the pion decay constant /,r,ph — 92.4 MeV 
are fixed to their empirical values, since they refer to the 
parameterization of the nuclear (matter) background. On 
the other hand, the parameter F^^ — 2/7r of the Skyrme 
Lagrangian is fixed to the value 108.11 MeV, since this pa- 
rameter together with e = 4.826 and m.^ = 134.977 MeV 
refers to the soliton itself; i.e. the empirical (isospin-ave- 
raged) masses of the nucleon and delta are reproduced by 
this choice, as explained below Eq. ([4]). The value of the 



cutoff-scale parameter Ac = 737 MeV is adjusted to the 
the threshold value of the isospin-odd on-shell ttN scatter- 
ing amplitude. The dimension-two low-energy constants 
Ci,2,3 can be found e.g. in Refs. [37ll38]. These values are 
consistent with the most recent pion-nuclei scattering data 
as it was summarized in Ref. [47] and lead to the thresh- 
old values 6o("^7r) ~ and bi{m.,r) « — 0.0883to~^, respec- 
tively . Within the errors these values of the scattering 
lengths are consistent with the more refined analysis of 
Ref. [49]. Furthermore, the incorporation of double scat- 
tering corrections in the s-wave pion polarization operator 
leads to the effective isoscalar scattering length 



.off 



(a;)«6o(^)-^[6g(^) + 26?(a;)], (22) 



where kp = [37r2p/2]^/'^ is the total Fermi momentum. 
The terms of higher order than cj^ can be neglected in 
and 61, provided that the condition uj < m-,^ is met(3 

For simplicity, we ignore the energy dependence in the 
p-wave scattering volumes and replace Co(w) and Ci(a;) 
by the constant threshold values co(TO7r) = 0.21to~3 and 
Ci(m7r) = 0.165m^'^ of the 'current' SAID analysis |50j . 
This is compatible with the discussion in Ref. [41]. Fur- 
thermore, all terms proportional to odd uj powers in 77° 

and even ones in AIIp are neglected. This is consistent 
with the remark in footnote [D and the disregard of pion- 
absorption in this approach. 

In summary, one can write the polarization operators 
([TU)1 as follows: 



71° = 



Xs {^)+xt{^) 



~ 'ikp ~9 
^0 + TT^K P 



bo 



3kp 



02 -2 2 

Xs ^ ^ 



11'' 

-L-Lp 



= Vx° • V « V 47rco(TO^)p/7y 

^ 1 + 4Trg'co{m^)p/ri 

Alls = -bi5pm:j^^LU = -Axs'm'^^^uj , 



V, 



An„ 



27: UJ 
mNT] 



ci(m^) (V^Sp) = -Axpvn^^uj . 



(23) 



(24) 

(25) 
(26) 



2.3 Medium-modified Lagrangian 

Evidently the explicit expressions of the self energies in 
configuration space follow from the standard rules 77° (tj) 

-> 77°(i9o) and An{uj) An{ido). After inserting the 
polarization operators n{ido) and AII{ido) from ([23)) - ([26)) 
into the Lagrangian (111!) and integrating by part, in order 



^ Within the framework of the Skyrme model, this situation 
corresponds to nucleons with S — T — ^ '-^ loA, where A ~ 
1 fm is the moment of inertia of the skyrmion. In the case of 
Zi-isobar states (5 = T = | ~ uA) also uj" terms with n > 3 
have to be taken into account. 
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to symmetrize in the time derivatives, one arrives at the 
final form of the medium-modified Lagrangian: 



sym 

r* — AC 



XSB 



(27) 
(28) 
(29) 



16 



-(l-X°)Tr(VC/.V[/t)}, (30) 

p2 2 

^XSB - 16"*^ ' 

xTr[(C/-l)(C/t-l)] , (31) 
AC^.s = -^^Al'_Tr(r„C/)Tr(raC/^), (32) 

a— 1 



a,b=l 



X Tr(T,C/)Tr(Tb9ot/^) ■ 



(33) 



Here AC^aes and ^vC|nv ^''^^ the isospin-breaking terms 
arising from the explicit symmetry breaking in the mesonic 
sector and the isospin asymmetry of the surrounding en- 
vironment, respectively. 

Note that both the temporal part of £2 and the chi- 
ral symmetry breaking term 'C*gg, decrease - at lead- 
ing order linearly - with increasing matter density, since 
and ^'''^ negative, see Eqs. (|20|) and ([23l) . How- 
ever, as the same equations indicate, ~ ~ ^oP, 
such that the effective mass, determined by the mass pole 
of the in-medium propagator, is approximately density- 
independent in agreement with the findings about the in- 
medium Gell-Mann-Oakes-Renner relation of Refs. |45l 
fiS] , Furthermore, one can see that the Lagrangian (|33p 
contains the Weinberg- Tomozawa term, as the relation 



3 Classical solitonic solutions 

By following the two-stage method of Ref. [SSjfl of con- 
strained and unconstrained collective isospin-rotations, ap- 
plied to the hedgehog ansatz U = exp [ir • (r/r)F(r)], 
the time-dependent Lagrangian can be constructed from 
Eq. (|27p in terms of the standard angular velocities uoi of 
the collective modes and the constrained angular velocity 
a* (see below) as 



L* ^ I C*d^r = -Af*p - MiA^ + —A 



uj3(a*A* + A* 



a* i —A* + A' 



(35) 



Here 



A/, 



NP 



2 



252 



2F^ 



00 



+ Fl {ml + x°°) (1 - cos F) |r2 dr (36) 

is the in-medium mass of the soliton when it is not per- 
turbed (NP) by any isospin breaking. The abbreviations 
Fr = dF/dr and S = smF have been used, where F = 
F{r) is the chiral profile function of the hedgehog ansatz. 
Furthermore 



A* = {l + m-\f)A^+A, 

with the separate contributions 

00 

A_ = '^Fl I S^r^dr. 



F^ + ^]S^r^dr 



(37) 

(38) 
(39) 



b[°Sp 



(34) 



is based on the isovector s-wave scattering length in the 
chiral expansion to lowest order [5T|[52]. 

The Lagrangian (p7)l - p3)) will be used in our studies of 
isospin breaking effects in asymmetric nuclear matter. In 
the next sections we will present and discuss the changes 
that emerge due to the isospin asymmetry of the surround- 
ing nuclear environment. Specifically, we will concentrate 
on isospin-breaking effects in infinite nuclear matter with 
a constant density, so that the p-wave contribution pro- 
portional to Axp ^ 'V^Sp vanishes. Note that in the case 
of finite nuclei this term may be essential for nucleons lo- 
cated near the surface of the nucleus. 



is the in-medium moment-of-inertia, whereas 

A* = {2m,)-^AxsA^ (40) 

is the response of the isospin-asymmetric environment (see 
Eqs. (|25p and ((33)) ) to the collective iso-rotations. 

The constrained angular velocity parameter a* corre- 
sponds to a stationary rotation around the third axis in 
isotopic space that serves to undo the effect of the mesonic 
isospin-breaking term proportional to A4- at the classical 
level, when the collective rotational modes in the isospin- 
space are frozen (ct'1,2,3 0). In this classical limit, ap- 
plying the constraint [35j 



a 



*2 



2M'iA^/A* 



(41) 



An alternative, but equivalent 
method is presented in Appendix [A] 



way of introducing this 
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one generates the Lagrangian 

L* = -M^p + a* A* 



(42) 



The pertinent equation of motion for the hedgehog profile 
function F{r) takes then the form 

F^{l^xl){r^F.rr + 2rFr-S2) 



2S Fr\ 



S2 F, 



52 



Fl{ml + xT)Sr^ + 



3m^ 



527-2 = 0, (43) 



where the additional abbreviations S2 ~ sin 2F and Frr 
= (PF/dr^ were introduced. 

The solution corresponding to the soliton of baryon 
number B = 1 fulfills the boundary conditions 



lim F{r) = tt — Cr , 

lim F{r) — D {1 + mpr) exp {— m^r} /r^ , 

r — ^CJO 

2a*m-Mxs/3 



ml 



ml + xT 



(44) 
(45) 

(46) 



where C and D are constants. Since the parameter a* is 
part of the classical equation (|43|) . i.e. h{Frr,Fr,F,a*) = 
0, this equation together with the constraint (|4ip can be 
solved by iteration: 



h{F^?,F^''\F(°\ 0)^0 



h{F^:^> , Fr , , a; _i) = => < = a* {F^ , . 

In the actual calculation, this iteration scheme rapidly 
converges after 3 to 4 iteration steps. 



4 In-medium neutron-proton mass difference 

4.1 Strong part of Am*^ 



By taking into account the condition (j4T|) . applying the 
definition a*ff = a* + A* / A* and using the canonical 
quantization procedure as in Ref. [35 , one can construct 
from the Lagrangian (I35|) the Hamiltonian in terms of the 
isospin operator T: 



H* = AF 



NP 



T 



Ti T2 

2A* 2A* 
2 



T3 - A*a; 



= A/^p + ali^n + A* ^ 



(47) 



Thus the strong part of the neutron-proton mass difference 
can be identified as 



A* 

A* 



(48) 



The choice of the sign of a* is fixed by the sign of Zim^p""^^ 
in free space [35] . 



Note that the density-variation of the strong part of the 
neutron-proton mass difference will be more pronounced 
than in Ref. [35 for the following reasons: (i) the ex- 
plicit density-dependence of the moment of inertia A* (see 
Eq. (|37p) resulting from the energy-dependent parame- 
terization of the optical potentials, and (ii) the existence 
of the additional term A* / A* in an isospin- asymmetric 
environment. Even if the explicit isospin breaking in the 
mesonic sector were omitted, = 0, there still would 
be a non- vanishing neutron-proton mass splitting propor- 
tional to the isospin-asymmetric environment factor A* . 



4.2 Electromagnetic part of Am^^ 

As discussed in Appendix [b1 by calculating the pertinent 
Noether currents one can construct the following isoscalar 
(S) and isovector (V) electromagnetic (EM) form factors 



B A* \ 
--—A + A*\joiqr) dr. 



m^il + A*) f~^2Jiiir) 



2A* 



qr 



Gr(g2) 

Gl*{q')=m^ I (l-xl)A^+Al 



1 

2A* 



qr 



(49) 
(50) 
(51) 
dr , (52) 



in terms of the spherical Bessel functions jo and ji and the 
three-momentum transfer q = \q\. Here a quantity with a 
tilde, say Z = Z{F), is defined as the integrand of the 
corresponding functional, i.e.: 

00 

Z[F] EE J Z{F{r)) dr. 



As usual, _B = 1 is the baryon charge, such that B{r) 
= 4TTr'^B°{r), where B°(r) = - sin^ F^;/(27r2r2) is the 
baryon density of the skyrmion. The medium-dependent 
form factors of the proton and neutron are defined as 



g|:m(9^)±g^:m('7' 



with the normalization conditions Gg*(0) ~ 1, Gg*(0) = 
0, G^(0) = ^p, Gf^{0) = /X*, where /i* and /i* are the 
magnetic moments of the in-medium proton and neutron, 
respectively. 

In the present approach all form factors explicitly de- 
pend on medium functionals, on the one hand, via the 
density-dependent moment of inertia A* (see Eq. ([57)) ). 
and on the other hand, via additional terms resulting 
from the isospin-asymmetric nuclear environment. More- 
over, note the additional terms in the isoscalar form fac- 
tors as compared with Ref. |35| . which emerge here from 
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that part of the isospin charge density that is independent 
of the isospin T3 (see Appendix IB|) . 
Finally, applying the formula |53| 



2m^ 



(53) 



where a « 1/137 is the fine-structure constant, one can 
calculate the medium-dependent electromagnetic part of 
the neutron-proton mass difference as in Ref . |35j . 




0.4 0.6 
P/Po 



5 Results and discussions 

In Fig. [1] the strong part of the in- medium neutron- proton 
mass splitting, z4m*p, is shown for isospin-symmetric nu- 
clear matter (solid curve), neutron-rich matter (dashed 
curve), pure neutron matter (dotted curve), and proton- 
rich matter (dot-dashed curve) . Already in isospin-symme- 



> 

(U 



o 




0.4 0.6 
P/Po 



Fig. 1. Density dependence of the strong part Am 



.(strong) 



of 



the neutron-proton mass difference. The abscissa represents 
the density p normalized to the saturation density of ordinary 
nuclear matter po = 0.5m^, while the ordinate shows the mass 
difference in units of MeV. The result in isospin-symmetric 
matter is plotted as a solid curve, the result of neutron-rich 
matter with 5p/p — 0.2 as dashed curve, the dotted curve 
represents pure neutron matter {5p/p = 1) and the dot-dashed 
curve proton-rich matter with 6p/ p = —0.2. 



trie matter At 



* (strong) 



has visibly a different density- 



behavior than the corresponding quantity of Ref. [35j . For 

example, at normal nuclear matter density, Z\mnp*'°'^^^ 
has increased by about 42% relative to its free space value 
(see the solid curve of Fig. [1]). In contrast to this in the 
previous work |35j . where the optical potentials were as- 



sumed to be energy- independent, Ai 



* (strong) 



decreased 



Fig. 2. Density dependence of the electromagnetic part 
Am'^\^^^ of the neutron-proton mass difference. The axes and 
curves are defined as in Fig. [T] 



nuclear matter densities; in other words, Zlrrinp'™"®^ in 
Ref. 35J was practically density- independent. 

Moreover, when the isospin symmetry of nuclear mat- 



ter is broken, Sp/p ^ 00 Ai 



* (strong) 



strongly varies at the 



qualitative as well as quantitative level (see the dashed 
curve in Fig. [1]). In pure neutron matter the change be- 
comes very drastic (see the dotted curve in Fig. [Ij, and 
^?Tii*p*'°'^^^ decreases very fast with increasing density. 

In contrast to the strong part, the electromagnetic part 
of the neutron-proton mass difference varies only by a 
small amount when the isospin-asymmetry parameter is 
increased (see Fig. But compared with the result of 
the previous approach I35J, the change is sizable, even in 
isospin-symmetric matter. This is again due to the explicit 
density dependence of the moment of inertia (|37[) , and the 
changes in the solutions of the classical equations 
Note that with increasing density the moment of inertia 
A* decreases since Xs^ < 0- In addition, the solutions of 
the classical equations are altered because Xs > in 
Eq. (22) of Ref. [3S] is replaced by the w-independent part 
of the present Eq. ([23)1 . namely by which is negative. 

For completeness, we present the total neutron-proton 
mass difference in Fig. [31 From a comparison with Fig. [T] 
it is obvious that this mass difference is completely dom- 
inated by its strong part. In pure neutron matter and 
at the density pq, the neutron-proton mass difference is 



Ami 



-25 MeV. For comparison, the authors of the 



work [IT] got the result ^m*p — 70 MeV in framework 
of QCD sum rule studies. 

Another interesting result is the difference between the 
values of ^'7i*p in neutron-rich and proton-rich matter - 
compare the dashed and dash-dotted curves in Fig. [31 One 
can see that in neutron-rich matter Am^p is decreased rel- 
ative to the isospin-symmetric case, whereas in proton-rich 
matter the behavior is just opposite. This finding may be- 
come useful for future studies of mirror nuclei and their 



by a very tiny amount, namely by about 2% at normal 



^ This quantity may be called the isospin-asymmetry param- 
eter of the nuclear environment. 
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Fig. 3. Density dependence of the total neutron-proton mass 
difference /im*p. The axes and curves are defined as in Fig. [T] 



properties. For example, for the case of the mirror nu- 
clei ""^Ca and ''^Ni, a similar behavior of /im*p was found 
in Ref. 30J within density-dependent relativistic hadron 
field theory. The result of our work is also consistent with 
the findings of Refs. |28l[24l[25l[26l[27l[28l[29] that utilize a 
relativistic approach and with the nonrelativistic calcula- 
tion [21] based on Skyrme-like effective interactions. 

Even at the qualitative level, the various models men- 
tioned in the introduction differ in their predictions of 
the neutron-proton mass difference in nuclear matter: (i) 
in nonrelativistic approaches |18, 19,20J, which are focused 
on the system properties as a whole, this difference mainly 
turns out to be positive (Z\m*p > 0); (ii) however, it is 
negative (Z\to*p < 0) in relativistic approaches p3 l[24ll25l 
[26l[27[|28l[29] and some nonrelativistic variational calcula- 
tions [H] or it becomes negative with increasing density 
in QCD sum rule studies [^ fTUlfTT] : (iii) it depends on the 
isospin content of the system (/4m*p > or /\m*p < 0) in 
relativistic hadron field theory [30] . The effective masses in 
relativistic approaches are discussed in detail in Ref. [S3^ . 
Furthermore, the difference in the behavior of ^?7i,*p in the 
relativistic and nonrelativistic approaches is explained in 
Ref. [H]. 

Also our approach shows that Zim*p qualitatively de- 
pends on the isospin content of surrounding environment. 
It is always positive in proton-rich matter as well as in 
isospin-symmetric matter (see the solid and dot-dashed 
curves in Fig. [3]). In neutron- rich matter, however, the 
sign may change. For the reader's convenience, we plot 
in Fig. d] those values of the density p as function of the 
isospin-asymmetry parameter Sp/p where the in-medium 
neutron-proton mass splitting or its strong part vanishes. 
In other words the solutions of the 2-parameter equations, 
Z\m*p(p, = and Amnp^™^\p,5p) = 0, are pre- 
sented. For small positive input for the isospin-asymmetry 
parameter the neutron-proton mass difference or its strong 
part vanishes at high densities (of the order of the ordi- 
nary nuclear matter density po). With increasing Sp/p, 
however, the mass difference changes its sign at moder- 




0.2 0.4 0.6 0.8 1 



Sp/p 

Fig. 4. The solutions of the 2-parameter equations 
Z\m*p(/9, — and zim^p (p, 5p) = 0. The abscissa 
represents the isospin-asymmetry parameter, while the ordi- 
nate shows the density (in units of the ordinary nuclear matter 
density po = 0.5m^), where the neutron-proton mass difference 
(solid curve) or its strong part (dashed curve) vanishes. 

ate densities, and in strongly isospin-asymmetric matter 
this transition is already at low densities. For instance, in 
neutron-rich matter with the isospin-asymmetry parame- 
ter Sp/p ~ 0.1 the proton becomes heavier at the density 
p ~ G.85po- In pure neutron matter this change happens 
already at the density p ~ O.OQpo- 

In addition, in Tabled] we present the calculated effec- 
tive masses and isoscalar as well as isovector charge radii 
of the in-medium nucleons for some values of the nuclear 
matter density. |fj In general, the nucleon masses strongly 
decrease in the nuclear medium and are qualitatively in 
agreement with the well known results |55p56j. At nor- 
mal nuclear matter density and for an isospin asymmetry 
Sp/p ^ 0.25, the difference in the effective masses (normal- 
ized to the corresponding free space values) of the neutron 
and proton, respectively, is m*/mn — ra*^/m-p ~ 0.01. For 
comparison, the result of Ref. [3D] for nucleons located 
near the center of "'^^^Sn is one order of magnitude bigger: 

TO*/mn — TOp/rTlp ~ 0.1. 

The isoscalar and isovector charge radii {r'^)*^{^Q i in- 
crease with increasing density of the mediumj^ The iso- 
scalar electric radius is more strongly affected by the iso- 
spin asymmetric environment than the isovector one be- 
cause of the presence of the isospin breaking term A* (see 
appendix [B]) . Consequently, in proton-rich matter the ra- 
tio (?'^)e^i^o/('"^)i/i=o iTiore enhanced than in neutron- 
rich matter. 

^ Note that the tabulated values of the free proton and neu- 
tron mass differ from their PDG values [43] since the customary 
Skyrme value Mjv = 938 MeV was used here and in Ref. [35] 
as input for the isospin-averaged nucleon mass. 

^ Note that our results in free space differ from the ones 
of Ref. [57] by a factor \/2 due to different normalizations of 
the charge densities (see Eq. (|B.3P in the Appendix |B] and the 
corresponding definitions in Ref. [57]). 
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Table 1. Calculated masses (in units of MeV) and isoscalar as 
well as isovector charge radii (in units of fm) of the nucleons in 
nuclear matter of density p (in units of the saturation density 
of ordinary nuclear matter po = 0.5m^). 



p/po 


» » ; 2\*l/2 


, 2X.1/2 
\' /E,I^1 





In free space 
937.4 938.7 0.49 


0.74 


0.5 
1.0 


In proton-rich matter {Sp/p = 
729.6 733.3 0.61 
547.9 555.9 0.79 


-0.2) 
0.84 
0.98 


In isospin symmetric matter {5p/ p — 0) 
0.5 729.7 731.4 0.60 0.84 
1.0 547.9 550.2 0.75 0.98 


0.5 
1.0 


In neutron-rich matter (Sp/p = 
731.5 731.3 0.58 
553.4 550.0 0.72 


-- 0.2) 
0.84 
0.98 


0.5 
1.0 


In pure neutron matter (Sp/p 
757.9 750.1 0.54 
632.1 607.0 0.52 


= 1) 
0.83 
0.94 




Fig. 

The 
The 



0.4 0.6 
P/Po 

5. Density dependence of the proton magnetic moment, 
ordinate represents p* in nuclear Bohr magnetons (n.m.). 
other definitions are as in Fig. [1] 



The density-dependence of the magnetic moments of 
the in-medium proton and the neutron is presented in 
Fig. [5] and Fig.[6l respectively. The influence of the isospin 
asymmetry of the surrounding environment on the in- 
medium magnetic moments is comparatively weak in this 
case. Relative to the result in isospin-symmetric matter 
both, the proton and neutron magnetic moments are de- 
creased in neutron-rich matter and increased in proton- 
rich matter. 

Let us conclude by remarking that within the present 
approach the change of -^mj^p is completely dictated by 
its strong part when the isospin-asymmetry parameter is 
sizable (compare Figs. [U [2] and [3]). 




0.4 0.6 
P/Po 

Fig. 6. Density dependence of the neutron magnetic moment. 
The ordinate represents pi in nuclear Bohr magnetons (n.m.). 
The other definitions are as in Fig. [1] 



6 Summary and outlook 

We have proposed an effective Lagrangian which incorpo- 
rates energy-dependent optical potentials for the s- and 
p-waves and which takes into account the influence of the 
isospin-asymmetry of the environment onto the Skyrme- 
solitons of the model. As a result the neutron-proton mass 
splitting in asymmetric nuclear matter is predicted to vary 
strongly relative to its free space value. The predictions 
obtained in the present work are in a qualitative agree- 
ment with the ones obtained within relativistic hadron 
fleld theory [28[|29[[30] and nonrelativistic variational cal- 
culations [21]. Quantitatively, however, the changes of Z\m*p 
are small in comparison to the results of those works. Our 
approach shows that ^ml^^ in nuclear matter with siz- 
able asymmetry is mainly dictated by its strong part. In 
the case of more complicated calculations involving finite 
nuclei this may serve as a justification to evaluate only 
the strong part of the mass difference within the present 
approach. 

Consequently, the next step in our future studies will 
be the estimate of ^w^p in finite (particularly in mirror) 
nuclei. Here additional effects are expected because the 
p-wave contribution proportional to ~ ^^Sp in the 
Lagrangian (|33p and the deformation effects discussed in 
Refs. [58u59] become relevant. 

The work of U.T.Y. was supported by the Alexander von Hum- 
boldt Foundation. Partial financial support from the EU Inte- 
grated Infrastructure Initiative Hadron Physics Project (con- 
tract number RII3-CT-2004-506078), by the DFG (TR 16, 
"Subnuclear Structure of Matter") and by BMBF (research 
grant 06BN411) is gratefully acknowledged. 



A Two-stage method and classical solutions 

In this appendix we review a new interpretation for the 
inclusion of the stationary a* rotations that is different 
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from the one of Ref. [36] or the one of Ref. [35] , but which 
leads - under the same input - to the same results as in 
those references. This new strategy is more convenient for 
our analytic calculations. 

First of all we remark that stationary a* rotations es- 
sentially have to be introduced because of the explicit 
isospin breaking in the mesonic sector, i.e. when A4- ^ 0. 
If 7V1_ — then a* — 0. Of course, the condition (|iT|) 
satisfies this requirement. Secondly, at the classical level, 
at which the soliton solution is determined, this mesonic 
mass splitting can be ignored - at the quantum level, how- 
ever, this can not be the case since the symmetry break- 
ing effect is enhanced by the coupling to the collective 
zero modes. Note the quadratic dependence on the small 
parameter at the classical level, whereas at the quan- 
tum level the dependence is linear; compare Eq. (|4ip with 
Eq. dUl). 

Let us, for the moment, put C*^ = in the Lagrangian 
(|?7|) . Then one can construct the classical 'hedgehog' so- 

lution(s) from the static Lagrangian £sym'^*"^^. In the work 
[35] these solutions are referred to as the solutions of the 
non-perturbed (NP) system, determined by 



5F 



(A.l) 



where F{r) is the usual chiral profile function of the hedge- 
hog ansatz and where M^p = — / £sym^*"^^ d^r (compare 
with Eq. p5p ) is the in- medium mass of the soliton when it 
is not perturbed by any isospin breaking. The Lagrangian 



*(static) 



tat ions 



sym is invariant under time-independent isospin ro- 



U AUA^ , A = exp{iT ■ ip' /2} , 



(A.2) 



which can be treated as the usual zero modes of the model. 
The resurrection of the time dependence ip^ — > <f^{t) in 
Eq. (|A.2p leads to the (spatially integrated) Lagrangian 



^sym d'^r — — Afj^p 



A* 



where A* is the in-medium moment-of-inertia. 

Now we plug the isospin breaking 7^ back into 
the mesonic sector, such that the corresponding Lagrangian 
reads 



out when the angular integration is performed, although 
Ai- ^ 0. But in order to be consistent in the use of the un- 
changed classical solution F{r) determined by Eq. (|A.1[) . 
one has to introduce - at the classical level - the require- 
ment ^ 

I^^A* ~ MlA^ ^0. (A.5) 

One can see explicitly that classically not all the rotations 
Lpl,(P2, ^3 of Eq. (|A.2[) can be time-independent. In other 
words, at least one of the angular velocities (pl,ip2,^l 
must not vanish, say, the one with respect to the third 
axis in isotopic space (which will later be the quantiza- 
tion axis). We will call the corresponding (via Eq. (|A.5|1 ) 
constrained angular velocity a*, such that the condition 
(|A.5[) is nothing else than the constraint (|4T]) . As men- 
tioned above, after angular averaging, the system remains 
invariant under time-independent zero modes which we 
called (fl^ , ip2^ , ipl^ . Again these zero modes can be made 
time-dependent 



d^. 



i C* 

\ sy] 



= -M, 



/^2 



NP 



■A*-MlA^ 



(A.3) 

Note that the integrated Lagrangian (IA.3p . even in the 
presence M- ^ 0, is still invariant under additional time- 
independent isospin rotations, e.g. 

Ait) ^ A{t) = exp{ix • [ip\t) + ip'']/2} 

EE exp{ix • <p(t)/2} . (A.4) 

In fact, it is a property of the hedgehog ansatz that these 
time-independent isospin rotations can be compensated 
by time-independent spatial rotations. The latter average 



^1, 



^2, 



^3, 



(A.6) 



where lui,lu2,i-U3, in contrast to a*, are unconstrained an- 
gular velocities. Transcribed to our starting point (|A.4p . 
we eventually arrive at 



953 = ^3 + a* 



(A.7) 



Inserting ipi,ip2,(p3 into Eq. (|A.3P one generates the La- 
grangian ([55]) with A* = 0. Interpreting now 0i = (^i, (/)2 = 
f2,4'3 = f3—o.*t as collective coordinates and lui,lu2, lu^ as 
their pertinent velocities, the standard quantization pro- 
cedure should be performed around the stationary point 
a*t with the constant angular velocity a*. 

This way of handling the classical a* rotations is equiv- 
alent to the procedure performed in Refs. 35,36 . But 
there is no need anymore for the explicit introduction of 
the matrix T{t) (see equation (15) in Ref. [35]) and the 
interpretation of the corresponding matrix order ambigu- 
ities discussed in Ref. [36] . 

Thus at the analytical level one performs the usual 
Skyrme model calculations, where Eq. (|A.7p is inserted 
into the expression (IA.3P and where the constraint (|4T]) 
is applied at the classical level, such that L* = ~M^p 
is extremized under the hedgehog ansatz, see Eq. (jA.ip . 
The latter corresponds to classical equation of motion (|43]) 
for the iso-symmetric matter case with vanishing dp and 
therefore, according to the definition ([25]) , vanishing Axs 
(and vanishing A* , see Eq. ([40]) '). 

Finally, the inclusion of isospin-breaking effects due to 
the asymmetric environment {AC*^^^ 7^ 0) leads to the 
Lagrangian 



dV \ + AC 



NP 



-A* 



+ Ac: 

MlA^+ipA*. (A.8) 



Note that the Lagrangian (|A.8|) . because of the angu- 
lar averaging, remains invariant under time-independent 
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isospin-rotations as its counterpart (jA.Sp . The collective 
coordinates are therefore 0i ~ tpi, (f)2 = ^2,4''i = V5 ~ 
as before and Eq. (jA.7[) is still valid. When the latter is in- 
serted, Eq. (jA.8|l transforms to the Lagrangian ((35|) which 
is the starting point for the quantization in this work. 
But in order to be consistent with the use of a solution of 
hedgehog type, one still has to satisfy two requirements: (i) 
to apply the constraint ([1T|) in order to remove —Ai^A^ 
and A*{(p)'^/2 from (IA.8|) . and (ii) to extremize the re- 
mainder —M^p + (pA* (compare with Eq. (|42|) 1 under the 
hedgehog ansatz. Here (p{t) still corresponds to (p^t) = a* , 
i.e. if-^^ is still assumed to be time- independent. Conse- 
quently, instead of Eq. (lA.ip it is Eq. (gSl) with Axs i^ 
that has to be solved (in practice by iteration) in order to 
determine the profile function F(r). 



B Charges and magnetic moments 

In calculating the third component of the isospin current 
Vq'^\ one finds 



j dV^o^^) - (c^3 + a*) A* + A* = n 



(B.l) 



Consequently, the isospin charge density - modulo a factor 
47rr^ - is given as 

n = (W3 + a)A* +A* = (T3 -A*)^ + A*. (B.2) 

Note that there are terms that are independent of the 
isospin T3 on the r.h.s. Since the charges of the nucleons 
are defined as 

Q = ^+n = J pi=oir) d\ ± J pi^iir) dV , (B.3) 

the isoscalar and the isovector density distributions have 
here the following form 

Wp,^o(r) = I - ^A* + A* , (B.4) 



47rr^p/=i(r) 



2 A 
A* 
2A* 



(B.5) 



Analogous calculations for the magnetic moments lead to 



(P), ^ m^{l + A*) 



± 



6A* 



TTlN 



Br^dr 



{l~xl)A-+Al dr. (B.6) 
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